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The back-reaction of a classical gravitational field interacting with quantum matter fields is de- 
scribed by the semiclassical Einstein equation, which has the expectation value of the quantum 
matter fields stress tensor as a source. The semiclassical theory may be obtained from the quantum 
field theory of gravity interacting with N matter fields in the large N limit. This theory breaks 
down when the fields quantum fluctuations are important. Stochastic gravity goes beyond the semi- 
classical limit and allows for a systematic and self-consistent description of the metric fluctuations 
induced by these quantum fluctuations. The correlation functions of the metric fluctuations ob- 
tained in stochastic gravity reproduce the correlation functions in the quantum theory to leading 
order in an f/N expansion. Two main applications of stochastic gravity are discussed. The first, 
\Q , in cosmology, to obtain the spectrum of primordial metric perturbations induced by the inflaton 

f"^ . fluctuations, even beyond the linear approximation. The second, in black hole physics, to study the 

f^*) ' fluctuations of the horizon of an evaporating black hole. 
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I. INTRODUCTION 

t-H . 

There are two aspects to semiclassical gravity. On the one hand we have quantum field theory in a curved spacetime 
which is now a well understood and well defined theory both for free fields [l[ and interacting fields Q . In this theory 
the gravitational field is the classical field of general relativity, the metric of the spacetime, and the quantum fields 
are test fields which propagate in such a spacetime. Since the spacetime is now dynamical it is not always possible to 
define a physically meaningful vacuum state for the quantum field and when this is possible in some "initial" times it 
\Q \ is usually unstable, in the sense that it may differ from the vacuum state at latter times, and spontaneous creation of 
■ particles occurs. Applications of this in cosmology, such as particle production in expanding Friedmann-Robertson- 
""q 1 Walker models (HQ, and black hole physics, such as Hawking radiation @, 0], are well known. This is one aspect of 
^jH, the interaction of gravity with quantum matter fields 

The second aspect of this interaction is the back-reaction of the quantum fields on the spacetime. Since the 
(3JT)! gravitational field couples to the stress tensor of matter fields, the key object here is the expectation value in a 
given quantum state of the stress tensor of the quantum field, which is a classical observable. However, since this 
object is quadratic in the field operator, which is only well defined as a distribution on spacetime, it involves ill 
defined quantities which translate into ultraviolet divergences. To be able to define a physically meaningful quantity a 
regularization and a renormalization procedure is required. The ultraviolet divergences associated to the expectation 
value of the stress tensor are also present in Minkowski spacetime, but in a curved background the renormalization 
procedure is more sophisticated since it needs to preserve general covariance. A regularization procedure which is 
specially adapted to the curved background is the so called point-splitting method @, H, Q . The final expectation 
value of the stress tensor using point splitting or any other reasonable regularization technique such as dimensional 
regularization, is essentially unique, modulo some terms which depend on the spacetime curvature and which are 
independent of the quantum state. This uniqueness was proved by Wald [ToL [Til Il2| who investigated the criteria 
that a physically meaningful expectation value of the stress tensor ought to satisfy. 

The back- reaction is formulated in terms of the semiclassical Einstein equations. These are Einstein equations 
which have the expectation value in some quantum state of the stress tensor as a matter source. The back-reaction 
problem was investigated in cosmology, in particular to see whether cosmological anisotropics could be damped by 
back-reaction [TH, [l4j • This was an earlier attempt [l5|, [l6[ previous to inflation to explain why the universe is so 
isotropic at present. 

A useful approach to the back-reaction problem is to use the effective action methods [13, EH Gil that are familiar in 
quantum field theory. These methods were of great help in the study of cosmological anisotropies since they allowed 
the introduction of familiar perturbative treatments into the subject. The most common effective action method, 
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however, led to equations of motion which were not real because they were tailored to compute transition elements of 
quantum operators rather than expectation values. Fortunately the appropriate technique had already been developed 
by Schwinger and Keldysh [2(| EI H^j in the so called Closed Time Path (CTP) or "in-in" effective action method, 
and was soon adapted to the gravitational context [23|, [24], [25], [2|| . These techniques were then applied to different 
problems in the cosmological context including the effects of arbitrary perturbations on homogeneous backgrounds 
[271 ] . As a result one was able to deduce the semiclassical Einstein equations by the CTP functional method: treating 
the matter fields as fully quantum fields and the gravitational field as a classical field. 

The semiclassical Einstein equations have limitations even outside the Planck scales: when the fluctuations around 
the exp ectation value of the stress tensor of the matter fields are large the semiclassical equations should break down 
[28l [23. [30j . One expects, in fact, that a better approximation would describe the gravitational field in a probabilistic 
way. In other words, that the semiclassical equations should be substituted by some Langevin-type equations with a 
stochastic source that describes the quantum fluctuations. A significant step in this direction was made by Hu [3lj who 
proposed to view the back-reaction problem in the light of the open quantum system paradigm, where the quantum 
fields play the role of the "environment" and the gravitational field plays the role of the "system" . Following this 
proposal a systematic study of the connection between semiclassical gravity and open quantum systems resulted in the 
development of a framework where different semiclassical Einstein-Langevin equations were derived [32l . l33l [34], |35| . 
The key technical factor to most of these results was the use of the influence functional method of Feynman and 
Vernon [36l . \3l\ for the description of the system-environment interaction when only the state of the system is of 
interest. The CTP method for open systems involves, in fact, the influence functional. 

However although several Einstein-Langevin equations were derived, these were always partial and related to some 
particular cosmological situation. On the other hand, the results were somewhat formal and some concern could be 
raised on the physical reality of the solutions of the stochastic equations for the gravitational field. This is related to 
the issue of the environment induced quantum to classical transition. In fact, for the existence of a semiclassical regime 
for the dynamics of the system one needs two requirements, in the language of the consistent histories formulation 
of quantum mechanics [38], H^]. The first is decoherence, which guarantees that probabilities can be consistently 
assigned to histories describing the evolution of the system, and the second is that these probabilities should peak 
near histories which correspond to solutions of classical equations of motion. The effect of the environment is crucial 
on the one hand to provide decoherence [40| and on the other hand to produce both dissipation and noise to the system 
through back-reaction, thus inducing a semiclassical stochastic dynamics on the system. As shown by Gell-Mann and 
Hartle [4lj in an open quantum system stochastic semiclassical equations are obtained after a coarse graining of the 
environmental degrees of freedom and a further coarse graining in the system variables. That this mechanism could 
also work for decoherence and classicalization of the metric field was not so clear lacking a full quantum description 
of the gravitational field, and the analogy could be made only formally [12] ■ 

An alternative axiomatic approach to the Einstein-Langevin equations which was independent of the open system 
analogy was introduced: it was based on the formulation of a general and consistent dynamical set of equations for a 
perturbative correction to semiclassical gravity able to account for the lowest order quantum stress tensor fluctuations 
of matter fields [43j. It was later shown that these same equations could be derived, in this general case, from the 
influence functional of Feynman and Vernon in which, the gravitational field is treated as a classical field and the 
quantum fields are quantized, the first being, in fact, the "system" and the seconds the "environment" [44[. Also, 
insp ired by results in some simple open quantum systems (45j and results of stochastic gravity in Minkowski spacetime 
[46l ] . the concern on the reality of the stochastic solutions was latter clarified. It was realized that the correlation 
functions of the metric fluctuations obtained in stochastic gravity reproduce the correlation functions in the quantum 
theory of gravity interacting with N quantum fields to leading order in an 1/N expansion [4?], [48| . Thus, stochastic 
gravity may be understood as a powerful and useful framework to study quantum metric fluctuations. 

Here we review the development of stochastic gravity and some of its applications. In section [IT] a brief sketch of 
semiclassical gravity is given. In section IIIII the axiomatic approach to stochastic gravity is discussed by introducing 
the Einstein-Langevin equations. In section IIVI to illustrate the relation between the semiclassical, stochastic and 
quantum theories, a simplified model of scalar gravity interacting with N scalar fields is considered. In section [V] 
we review an important application of stochastic gravity in cosmology. It concerns the computation of the two-point 
correlations of the metric perturbations induced by the fluctuation in the stress tensor of the inflaton field during 
inflation. The results to linear order agree with the standard results but the present method, in which the matter fields 
are treated exactly, may go beyond the usual approaches where the inflaton fluctuations are treated at linear level 
only. In section [VI] we deal with another important application in black hole physics: the study of the fluctuations 
near the horizon of an evaporating black hole. This subject is still under consideration and we only sketch some of 
the recent results. Finally, in section fVIII we summarize our result and briefly discuss other applications. We should 
mention that two reviews of stochastic gravity and its applications are now available [4^ . |50| . 
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II. SEMICLASSICAL GRAVITY 

Semiclassical gravity is a theory which describes the interaction of the gravitational field assumed to be a classical 
field with matter fields which are quantum. It is supposed to be some limit of the still unknown quantum theory 
describing the interaction of gravity with other fields. Due to the lack of the full quantum theory, the semiclassical 
limit cannot be rigorously derived. However, it can be formally derived in several ways. One of them is the leading- 
order 1 /N approximation of quantum gravity j5lj , where N is the number of independent free quantum fields which 
interact with gravity only. In this limit, after path integration one arrives at a theory in which the gravitational field 
can be treated as a c-numbcr (i.e. is quantized at tree level) and the quantum fields are fully quantized. If we call g a {, 
the metric tensor and <p the scalar field (for simplicity we consider just one scalar field) one arrives at the semiclassical 
Einstein equation as the dynamical equation for the metric g a b- 

G a b[g] = n(f ab [g]) ren , (1) 

where k = 8ttG = 8Tr/m'p, G is Newton's gravitational constant and rap is the Planck mass, T a b — T a b[4> 2 } is the stress 
tensor operator which is quadratic in the field operator <fr. This operator, being the product of distribution valued 
operators, is ill defined and needs to be regularized and renormalized, the subscript ren means that the operator has 
been renormalized. The angle brackets on the right hand side mean that the expectation value of the stress tensor 
operator is computed in some quantum state, say compatible with the geometry described by the metric g a b- On 
the left hand side G a b[g] stands for the Einstein tensor for the metric g a b together with the cosmological constant 
term and other terms quadratic in the curvature which are generally needed to renormalize the stress tensor operator. 
The quantum field operator <p propagates in the background defined by the metric g, it thus satisfies a Klein-Gordon 
equation, 

(V 2 g - m 2 - £%]) = 0, (2) 

where V 2 , stands for the D'Alambert operator in the background g a b, £ is a dimensionless coupling parameter (£ = 
is the minimal coupling and £ = 1/6 is the conformal coupling) and R is the Ricci scalar for the background metric. 
Equation |T]) is the semiclassical Einstein equation, it is the dynamical equation for the metric tensor g a b and describes 
the back-reaction of the quantum matter fields on the geometry. A solution of semiclassical gravity consists of the set 
(gab, 4>, \ip}) where g a b is a solution of ([I}, is a solution of © and is the quantum state in which the expectation 
value of the stress tensor in equation (fT|) is computed. 

This theory is in some sense unique as a theory where the gravitational field is classical. In fact, the (classical) 
gravitational field interacts with matter fields through the stress tensor, and the only reasonable c-number stress 
tensor that one may construct [ici| with the operator T ab is just the right hand side of (fT]), modulo the curvature terms 
needed for renormalization. However the scope and limits of the theory are not so well understood as a consequence 
of the lack of the full quantum theory. It is assumed that the semiclassical theory should break down at Planck scales, 
which is when simple order of magnitude estimates suggest that the quantum effects of gravity cannot be ignored: 
the gravitational energy of a quantum fluctuation of energy in a Planck size region, determined by the Heisenberg 
uncertainty principle, is of the same order of magnitude as the energy fluctuation. 

There is also another situation when the semiclassical theory should break down, namely, when the fluctuations 
of the stress tensor are large. This has been emphasized by Ford and collaborators. It is less clear how to quantify 
what a large fluctuation here means and some criteria have been proposed [30l l52l [53| . Generally this depends on the 
quantum state and may be illustrated by the example used in ref. [29( as follows. 

Let us assume a quantum state formed by an isolated system which consists of a superposition with equal amplitude 
of one configuration with mass Mi and another with mass M%. Semiclassical theory as described in (fTJ) predicts that 
the gravitational field of this system is produced by the average mass (M\ + M^fl, that is a test particle will move 
on the background spacetime produced by such a source. However one would expect that if we send a succession of 
test particles to probe the gravitational field of the above system, half of the time they would react to the field of a 
mass Mi and the other half to the field of a mass M2 . If the two masses differ substantially the two predictions are 
clearly different, note that the fluctuations in mass of the quantum state is of the order of (Mi — M2) 2 ■ Although 
the previous example is suggestive a word of caution should be said in order not to take it too literally. In fact, if 
the previous masses are macroscopic the quantum system decoheres very quickly [40( and instead of a pure quantum 
state it is described by a density matrix which diagonalizes in a certain pointer basis. For observables associated to 
this pointer basis the matrix density description is equivalent to that provided by a statistical ensemble. In any case, 
however, from the point of view of the test particles the predictions differ from that of the semiclassical theory. 
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III. STOCHASTIC GRAVITY 

The purpose of semiclassical stochastic gravity, or stochastic gravity for short, is to be able to deal with the situation 
of the previous example in which the predictions of the semiclassical theory may be too rough. Consequently, our 
first point is to characterize the quantum fluctuations of the stress tensor. 

The physical observable that measures these fluctuations is related to the two-point stress tensor correlations. Let 
us consider the tensor operator t ab = T ab — (T ab )I, where / is the identity operator, and introduce the noise kernel as 
the four index bi-tensor defined as the expectation value of the anticommutator of the operator t ab : 

N a bcd( x ^y) = 7j({tab(x),tcd(y)})- (3) 

This expectation value is taken in the background metric g a b which we assume to be a solution of the semiclassical 
equation (fT]). An important property of the symmetric bi-tensor, N abo d(x,y) = N c dab(y x) , is that it is finite because 
the tensor operator t ab is finite since the ultraviolet divergences of T ab are cancelled by the substraction of (T ab ). Since 
the operator T ab is selfadjoint N abcc i(x,y), which is the expectation value of an anticommutator, is real and positive 
semi-definite. This last property allows for the introduction of a classical Gaussian stochastic tensor ^ ab defined by 

(U(x)) s = 0, (U(x)U(y))s = N abcd (x, y). (4) 

This stochastic tensor is symmetric ^ ab = £ ba and divergenceless, V a £ a ;, = 0, as a consequence of the fact that 
the stress tensor operator is divergenceless. The subscript s means that the expectation value is just a classical 
stochastic average. Note that we assume that £ a b is Gaussian for simplicity, in order to include the main effect. The 
idea now is simple, we want to modify the semiclassical Einstein equation (fT]) by introducing a linear correction to 
the metric tensor g ab , such as g a b + h ab , which accounts consistently for the fluctuations of the stress tensor. The 
simplest equation is obtained by adding to the right hand side of equation (fT]), but written for the perturbed metric, 
the stochastic tensor just introduced. Substraction of the semiclassical equation (fTJ) leads to the following equation, 

G ( ab l9 + h} = [g + h]) ren + ^ ab [g}), (5) 

where we recall that the background metric g a b is assumed to be a solution of equation (fT]). As indicated by the 
superscript (1) this stochastic equation must be thought of as a linear equation for the metric perturbation h a b which 
will behave consequently as a stochastic field tensor. Note that the tensor t; a b[g] is not a dynamical source, since it has 
been defined in the background metric. Note also that this source is divergenceless with respect to the metric, and it 
is thus consistent to write it on the right hand side of the Einstein equation. This equation is gauge invariant with 
respect to diffeomorphisms defined by any field on the background spacetime [Hj|. If we take the statistical average, 
equation (O becomes just the semiclassical equation for the metric g a b + (h a b)s where the expectation value of T ab is 
taken in the perturbed spacetime. The quantum field now propagates in the spacetime described by the perturbed 
metric and thus it satisfies 

(V 2 g+h -m 2 - £R{g + h})4> = 0, (6) 

where the Ricci scalar is evaluated for the perturbed metric. 

The stochastic equation ([5]) is known as the Einstein-Langevin equation. The equation predicts that the gravitational 
field has stochastic fluctuations over the background g ab - It is linear in h ab , thus its solutions can be written as follows, 

h ab (x) = h° ab (x) + kJ dV ffiC^. x')C d (x'), (7) 

where h ab (x) is the solution of the homogeneous equation containing information on the initial conditions and 
G a \* cd (x,x') is the retarded propagator of equation (|5|) with vanishing initial conditions. Form this we obtain the 
two-point correlations for the metric perturbations: 

(h ab (x)h cd (y))s = (h° ab (x)h° cd (y))s + « 2 | d 4 x'd 4 yV5(^05(yOG^/(^^')^ e/9,l (2;',y')G^(y,2/')- (8) 

There are two different contributions to the two-point correlations. The first one is connected to the fluctuations of 
the initial state of the metric perturbations, we refer to them as intrinsic fluctuations. The second contribution is 
proportional to the noise kernel and is thus connected with the fluctuations of the quantum fields, we will refer to 
them as induced fluctuations. These two-point stochastic correlations are the most relevant physical observable, to find 
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them requires to know the noise kernel N a b c d(x,y). Note that the noise kernel should be thought of as a distribution 
function, the limit of coincidence points has meaning only in the sense of distributions. Explicit expressions of this 
kernel in terms of the two poin t Wightman functions is given in (43[, expressions based on point-splitting methods 
have also been given in [54l |55|. 

These stochastic correlations for the metric perturbations satisfy a very important property. In fact, it can be 
shown that they correspond exactly to the symmetrized two-point quantum metric correlations obtained in the large 
N expansion: 

^({hab(x),h cd (y)}) = (h ab (x)h cd (y)) s , (9) 

where h a b{x) mean the quantum operator corresponding to the metric perturbations. This result was implicitly 
obtained in the Minkowski background in ref. [46j where the two-point correlation in the stochastic context was 
computed for the linearized metric perturbations. This stochastic correlation exactly agrees with the symmetrized 
part of the graviton propagator computed by Tomboulis [56| in the quantum context of gravity interacting with N 
matter Fermion fields, where the graviton propagator is of order 1/N. This result can be extended to an arbitrary 
background in the context of the large N expansion [48| , a sketch of the proof with explicit details in the Minkowski 
background can be found in ref. [47J ■ This connection between the stochastic correlations and the quantum correlations 
was noted and studied in detail in the context of simpler open quantum systems [45] . It is thus clear that stochastic 
gravity goes beyond semiclassical gravity in the following sense. The semiclassical theory, which is based on the 
expectation value of the stress energy tensor, carries information on the field two-point correlations only, since (T a b) 
is quadratic in the field operator (j>. The stochastic theory on the other hand, is based on the noise kernel ([3]) which 
is quartic in the field operator. However, it does not carry information on the graviton-graviton interaction, which in 
the context of the large N expansion it gives diagrams of order 1/N 2 . This will be illustrated in section ITVl 



A. Functional approach 

To end this section we should mention that the Einstein-Langevin equation ([5]) may also be derived using the CTP 
functional method [13] • As remarked in the introduction the CTP functional was introduced by Schwinger [2(| [U, [22j 
to compute expectation values. One just considers the interaction of the gravitational field g a b, classical, with the 
field 4>, fully quantum. Then the effective action for the gravitational field is derived after integrating out the degrees 
of freedom of the quantum field, and the CTP influence action reduces basically to the Feynman and Vernon influence 
functional 36, used in quantum open systems. Here the system is the gravitational field and the environment 
is the quantum field. The stochastic terms for the gravitational field are found by suitably interpreting some pure 
imaginary term which appear in the influence action. These terms are closely connected to Gell-Mann and Hartle 
decoherence functional [4l| used to study decoherence and classicalization in open quantum systems. The net result 
is that the interaction with the environment induces fluctuations in the system dynamics. It is precisely the noise 
kernel introduced in ([3]) that accounts for this effect. 



IV. THE LARGE TV EXPANSION 



The large N expansion has been successfully used in quantum cromodynamics to compute some non-perturbative 
results. This expansion re-sums and rearranges Feynman perturbative series including self-energies. For gravity 
interacting with N matter fields it shows that graviton loops are of higher order than matter loops. To illustrate 
the large N expansion let us consider the following toy model of gravity, which we will simplify as a scalar field h, 
interacting with a scalar field <fr described by the Lagragian density 

L = \J d 4 x (d a hd a h + h(dh) 2 + ...)- J d A x (d a (f>d a (t> + m 2 <j) 2 ) + J d A x (h(dcj)) 2 + ...), (10) 

where, as previously k — 8ttG, and we have assumed that the interaction is linear in the (dimensionless) scalar 
gravitational field h and quadratic in the matter field (j) to simulate in a simplified way the coupling of the metric 
with the stress tensor of the matter fields. We have also included a self coupling graviton term of 0(h 3 ) which also 
appears in perturbative gravity beyond the linear approximation. 

We may now compute the dressed graviton propagator perturbatively as the following series of Feynman diagrams. 
The first diagram is just the free graviton propagator which is of O(k), as one can see from the kinetic term for the 
graviton in equation (|10|) . The next diagram is one loop of matter with two external legs which are the graviton 



6 



propagators. This diagram has two vertices with one graviton propagator and two matter field propagators. Since 
the vertices and the matter propagators contribute with 1 and each graviton propagator contributes with a n this 
diagram is of order 0(k 2 ). The next diagram contains two loops of matter and three gravitons, and consequently it is 
of order 0(k 3 ). There will also be terms with one graviton loop and two graviton propagators as external legs, with 
three graviton propagators at the two vertices due to the 0(h 3 ) term in the action (fT0)l . Since there are four graviton 
propagators which carry a k but two vertices which have n~ 2 this diagram is of order 0(n 2 ), like the term with one 
matter loop. 

Let us now consider the large N expansion. We assume that the gravitational field is coupled with a large number 
of identical fields <j>j, j = 1,...,N which couple only with h. Next we re-scale the gravitational coupling in such a 
way that R — kN is finite even when N goes to infinity. The Lagrangian density of this system is: 

N N 

L = — J d 4 x (d a hd a h + h{dh) 2 + ...)- J <?x (da^d"^ + m 2 cj) 2 ) + ^ J d 4 x {h{d^f + ...). (11) 

Now and expansion in powers of 1 /N of the dressed graviton propagator is given by the following series of Feynman 
diagrams. The first diagram is the free graviton propagator which is now of order 0(R/N) the following diagrams are 
N identical Feynman diagrams with one loop of matter and two graviton propagators as external legs, each diagram 
due to the two graviton propagators is of order 0(R 2 /N 2 ) but since there are N of them the sum can be represented 
by a single diagram with a loop of matter of weight N, and therefore this diagram is of order 0(R 2 /N). This means 
that it is of the same order as the first diagram in an expansion in l/N . Then there are diagrams with two loops of 
matter and three graviton propagators, as before we can assign a weight of N to each loop and taking into account the 
three graviton propagators this diagram is of order 0(R 3 /N), and so on. This means that to order l/N the dressed 
graviton propagator contains all the perturbative sums in powers of R of the matter loops. 

Next, there is a diagram with one graviton loop and two graviton legs. Let us count the order of this diagram: it 
contains four graviton propagators and two vertices, the propagators contribute as (R/N) 4 and the vertices as (N/R) 2 , 
thus this diagram is of 0(R 2 /N 2 ). Therefore graviton loop contributes to higher order in the l/N expansion than 
matter loops. Similarly there are N diagrams with one loop of matter with an internal graviton propagator and two 
external graviton legs. Thus we have three graviton propagators and since there are N of them, their sum is of order 
0(R 3 /N 2 ). To summarize, we have that when N — > oo there are no graviton propagators and gravity is classical, this 
is semiclassical gravity. To next to leading order, l/N, the graviton propagator includes all matter loop contributions, 
but no contributions from graviton loops and internal graviton propagators in matter loops. This is what stochastic 
gravity reproduces. 

That stochastic gravity is connected to the large N expansion can be seen from the stochastic correlations of linear 
metric perturbations on the Minkowski background computed in ref. (46j . These correlations are in exact agreement 
with the imaginary part of the graviton propagator found by Tomboulis in the large N expansion for the quantum 
theory of gravity interacting with N Fermion fields (56j . This has been proved in detail in ref. |47j and extended to 



the general case [48 1 



V. GRAVITATIONAL FLUCTUATIONS DURING INFLATION 



An important application of stochastic gravity is the derivation of the cosmological perturbations generated during 
inflation [54[ . Let us consider the Lagrangian density for an inflaton field <fi of mass m 

C^^^d^+^m 2 ^, (12) 

which is the basis of the simplest chaotic inflationary model [57| . The conditions for the existence of an inflationary 
period, which is characterized by an accelerated expansion of the spacetime, is that the value of the field averaged 
over a region with the typical size of the Hubble radius is higher than the Planck mass mp. After the Planck era one 
expects, by Heisenberg uncertainty principle, that in a region of the size of the Planck size nip 1 the energy density 
would be of the order of mp. If the inflaton field starts in one of those regions at a value much larger than mp then 
with a potential like that of eq. (|12[) one expects an inflaton mass m ^ mp in which case the potential is very flat 
and the field rolls slowly towards the minimum producing a period of almost de Sitter expansion of that region. In 
order to solve the cosmological horizon and flatness problem more than 60 e-folds of expansion are needed, to achieve 
this the scalar field should begin with a value higher than 3mp. Moreover, we will see that the large scale anisotropics 
measured [H,[5^, [6(| restrict the inflaton mass to be of the order of 10 -6 mp. 

We want to study the metric perturbations produced by the stress tensor fluctuations of the inflaton field on the 
homogeneous background of a flat Friedmann-Robertson- Walker model, described by the cosmological scale factor 
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a(rf), where rj is the conformal time, which is driven by the homogeneous inflaton field 4>(w) — {4>). Thus we write the 
inflaton field in the following form 

$ = 0(77) + <p(x), (13) 

where if(x) corresponds to a free massive quantum scalar field with zero expectation value on the homogeneous 
background metric: (ip) = 0. Restricting ourselves to scalar-type perturbations the perturbed metric g a b = gab + h a t 
can be written in the longitudinal gauge as, 

ds 2 = a 2 (?7)[-(l + 2$(x))dn 2 + (1 - 2^{x))8 i] dx l dx j ] 1 (14) 

where the metric perturbations ^>{x) and ^(x) correspond to Bardeen's gauge invariant variables [6lT [6^|. 

The Einstein-Langevin equation is gauge invariant, and thus we can work in a desired gauge and then extract the 
gauge invariant quantities. The Einstein-Langevin equation ([5]) reads now: 

G^[h}~K{f^[h})^ K U[9l (15) 

where g a b satisfies the semiclassical Einstein equations; assuming slow roll this background metric is an almost de 
Sitter metric. The superscript (1) refers to functions linear in the metric perturbation h a b- The stress tensor operator 
T a b for the minimally coupled inflaton field in the perturbed metric g a b is: 

fab = V Q 0V b + ^ a6 (V c 0V c + m 2 4> 2 ). (16) 

Now, using the decomposition of the scalar field into its homogeneous and inhomogcneous part, see eq. (|13[) . and 
the metric g a b into its homogeneous background g a b and its perturbation h a b, the renormalized expectation value for 
the stress tensor can be written as 

(fab[g}) = (f[g])<f«f, + (f ab [g])<f, v + {f ab [g}) vv , (17) 

where only the homogeneous solution for the scalar field contributes to the first term. The second term is proportional 
to (</?[<?]) which is not zero because the field dynamics is considered on the perturbed spacetime, i.e. this term includes 
the coupling of the field with h a b, as seen in equation ([6]). The last term corresponds to the expectation value to the 
stress tensor for a free scalar field on the spacetime of the perturbed metric. 

We can now compute the noise kernel N a bcd{x,y) defined in equation Q, which after using the previous decompo- 
sition may be written as 

({i a b,icd})[g} = ({i a bJcd})<j><p[g} + ({i a bJcd}) v <p[g}, (18) 

where we have used the fact that (0) = = (<p<p<p) for Gaussian states on the background geometry. We have 
considered the vacuum state to be the Bunch-Davies, or Euclidean, vacuum which is preferred in the de Sitter 
background, and this state is Gaussian. In the above equation the first term is quadratic in (p whereas the second one 
is quartic, both contributions to the noise kernel are separately conserved since both ^(77) and (p satisfy the Klein- 
Gordon field equations on the background spacetime. Consequently, the two terms can be considered separately. 
On the other hand if one treats ip as a small perturbation the second term in (jTSJ) is of lower order than the first 
and may be neglected consistently, this corresponds to neglecting the fluctuations associated with the last term in 
equation (|17[) . This approximation is equivalent to keep only linear terms in the inflaton perturbations. Stress tensor 
fluctuations due to a term like the last term of (jTTJ) were considered in ref. [631 ] . 

We can now write down the Einstein-Langevin equations (|15[) . It is easy to check that the space-space components 
coming from the stress tensor expectation value terms and the stochastic tensor are diagonal, i.e. (T,j) = = for 
i 7^ j. This, in turn, im plie s that the two functions characterizing the scalar metric perturbations are equal: $ = "J" 
in agreement with ref. 62]. The equation for $ can be obtained from the Oi-component of the Einstein-Langevin 
equation, which in Fourier space reads 

2ik i (H$ k + <f>' k ) = ^ F Z k0i , (19) 
mp 

where hi is the comoving momentum component associated to the comoving coordinate x l . Here primes denote 
derivatives with respect to the conformal time r\ and H = a! j a. A non-local term of dissipative character which comes 
from the second term in (fT7|) should also appear on the left hand side of equation (fT9")) . but we have ignored this 
term for simplicity (if one includes this non-local term it is then more convenient to write an equation combining the 
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other equations which is free of non-local terms; but the results are not substantially altered) . To solve this equation, 
whose left hand side comes from the linearized Einstein tensor for the perturbed metric [62|, we need the retarded 
propagator for the gravitational potential $> k , 



^V) = -^(%-,0^ + /(,,V)), (20) 



where / is a homogeneous solution of (|19[) related to the initial conditions chosen. For instance, if we take f(r], rf) = 
— ~ r l') a ( r l')/ a { r l) the solution would correspond to "turning on" the stochastic source at rjQ. 
The correlation function for the metric perturbations is now given by 



\2 



The correlation function for the stochastic source , which is connected to the stress tensor fluctuations through the 
noise kernel is given by, 

<& oMZ-k «(%)>, = \{{tUvi^?{m)})^ = ^M'M^'M^V,^), (22) 

where G k ^ (rji, 772) = ({ip k (r)i), <£>-fc(?72)}) is the fc-mode Hadamard function for a free minimally coupled scalar field 
which is in the Euclidean vacuum on the de Sitter background. 

It is useful to compute the Hadamard function for a massless field and consider a perturbative expansion 
in terms of the dimensionless parameter m/mp. Thus we consider g! (?71j V2) = a { r li) a (' r l2)G^ (rji, 772) = 
(Q\{yk(rii),y~k(r)2)}\0} = 27^(u fc (77i)Ufc(r? 2 )) with y k (rj) = a{rj)tp k {rj) = &kU k (j]) + <r_ k U*_ k (ri) and where u k = 
(2k)~ 1 / 2 e lkri (1 — i/rf) are the positive frequency fc-mode for a massless minimally coupled scalar field on a de Sitter 
background, which define the Euclidean vacuum state: a k \0) =0 Q. 

The background geometry, however, is not exactly that of de Sitter spacetime, for which a(if) = —(Hrj)^ 1 with 
— 00 < rj < 0. One can expand in terms of the "slow-roll" parameters and assume that to first order <j)(t) ~ m%(m/mp), 
where t is the physical time. The correlation function for the metric perturbation (f2Tj) is the computed, see ref. [54| 
for details. The final result, however, is very weakly dependent on the initial conditions as one may understand from 
the fact that the accelerated expansion of de quasi-de Sitter spacetime during inflation erases the information about 
the initial conditions. Thus one may take the initial time to be 770 = —00 and obtain to lowest order in m/mp the 
expression 

(**(r?)*k'(»?')>a ^Stt 2 (— ) k- 3 (27r) 3 S(k + k')co S k(7 1 -r ] '). (23) 
\m P J 

From this result two main conclusions are derived. First, the prediction of an almost Harrison-Zel'dovich scale- 
invariant spectrum for large scales, i.e. small values of k. Second, since the correlation function is of order of (m/mp) 2 
a severe bound to the mass m is imposed by the gravitational fluctuations derived from the small values of the Cosmic 
Microwave Background (CMB) anisotropies detected by COBE [58| and WMAP [59|,[6(|. This bound is of the order of 
(m/mp) ~ 10 -6 [58LI62I]. One possible advantage of the Einstein-Langevin approach to the gravitational fluctuations 
in inflaton over the approach based on the quantization of the linear perturbations of both the metric and the inflaton 
field [62j], is that an exact treatment of the inflaton quantum fluctuations is possible. On the other hand although 
the gravitational fluctuations are here assumed to be classical, the correlation functions obtained correspond to the 
expectation values of the quantum metric perturbations, as we have remarked in the previous section [4a. [48|. 

VI. FLUCTUATIONS NEAR BLACK HOLE HORIZONS 

Another interesting application of stochastic gravity is found in the context of black hole physics, in particular 
the stress tensor fluctuations near the black hole horizon may induce fluctuations in the horizon area. The relevance 



of these fluctuations in Hawking radiation needs to be understood 6J|. Some preliminary investigations seem to 
indicate that the fluctuations of the black hole horizon are always small and that the Hawking result should not be 
substantially different [H, Hi| , however, some other results by Bekenstein [67j seem to point in the opposite direction 
suggesting that the fluctuations of the black hole horizon may be significant in the long run. The contribution of the 
horizon fluctuations to the black hole entropy j68l , [69| is another interesting issue that may deserve some attention in 
the present context. 
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To clarify this situation Hu and Roura [7C| have analyzed this back-reaction problem in the stochastic gravity 
framework. Due to technical difficulties this is still work in progress so I will only summarize very briefly the main 
results. As shown by Hawking a black hole formed by spherical collapse emits thermal radiation with a temperature 
T = mp/8irM , where M is the mass of the black hole. This calculation was made under the assumption of quantum 
field theory in a curved background. That is, assuming that the black hole has a fixed mass much larger than the 
Planck mass mp, so that one can safely ignore quantum gravity effects, and that the black hole exterior can be 
described by the Schwarzschild metric. But, energy conservation arguments indicate that as the black hole emits 
radiation it will loss mass and evaporate. A precise calculation of the evaporation process requires the use of the 
semiclassical Einstein equations which describes the back-reaction of the quantum matter fields on the gravitational 
field in a self-consistent way. An exact self-consistent calculation of the evaporation process is by no means easy and 
has not been performed. Note that the black hole exterior is not vacuum any more, as the expectation value of the 
stress tensor (T a t,) ren is not zero, and is not described exactly by the Schwarzschild metric. Thus, even the radiation 
process needs to be reviewed. Furthermore even in the Schwarzschild background an exact analytic expression to 
describe that expectation value in the Unruh vacuum, which is the natural quantum state that describes the initial 
vacuum in a black hole formed by gravitational collapse, is not known. 

Fortunately, for large black holes the evaporation process is slow and a quasi-adiabatic approximation can be used to 
solve Einstein semiclassical equations. In this approximation one can assume that the black hole exterior is described 
by a Schwarzschild metric with a mass M which is the mass that the black hole has at that time. A suitable parameter 
to use in this approximation is the luminosity at a given time Lh = B/M 2 , where B is a constant which depends on 
the number of fields considered, the spins of these fields and the grey body factor. It has been estimated by Page [7l| 
to be of the order of 10~ 4 . Here, and in the rest of this section, we use units in which mp = 1. The quasi-adiabatic 
approximation holds as long as Lh <C 1. Black hole evaporation in the adiabatic approximation was described by 
Bardeen [72| and Massar [73| for spherically symmetric black holes. Let us summarize this calculation. A spherically 
symmetric metric can always be written as 



ds 2 = -e^ v ' r) 



^ 2m(v,r) 



dv 2 + 2e 4 ' (v > r) dvdr + r 2 (d6 2 + sin 2 0d<p 2 ). (24) 



There is an apparent horizon where the expansion of the outgoing null geodesies vanish, that is where dr a /dv = 0, 
which leads to r a (v) — 2m(v,r a (v)) = 2M(v). Thus we have defined 2M(v) as the apparent horizon, when M is 
constant this corresponds to the event horizon and M is the black hole mass. Thus in the adiabatic approximation 
we may consider that M(v) is the black hole mass at the "advanced" time v. The semiclassical Einstein equations 
become in the above coordinates, 

^WCO, ^ = -4«r 2 (T:), ?t = inr 2 {T rr ). (25) 
Ov Or ov 

We do not have an analytic expression for the expectation value of the stress tensor even in the Schwarzschild 
spacetime. However, at large radii it corresponds to a thermal flux of radiation and we may write (T£) = Lh /(Airr 2 ). 
Then one can use the stress tensor conservation equation to relate components on the horizon and far from it, 

W) +r2 gp = Q , (26) 
Or Ov 

Using this equation one may relate the positive energy flux radiated away from the horizon and the negative energy 
flux crossing the horizon. Taking this relation into account and the quasi-adiabatic approximation one finally gets, 
from the first of equations (|25p at the horizon, the equation for the evolution of the apparent horizon: 

dv ~ M 2 ' 10 

which gives the evaporation rate, as one would expect from energy conservation considerations. 

Now the Einstcin-Langevin equation may be used to study the metric fluctuations near the black hole horizon of 
the evaporating black hole. A full self-consistent computation is technically involved, in particular the computation of 
the noise kernel is very complicated, but Roura and Hu [7(| where able to give some reasonable estimates of the event 
horizon fluctuations in the evaporating process. They concentrate on the spherically symmetric fluctuations by pro- 
jecting the Einstein-Langevin equation on the spherical sector. The Einstein-Langevin equation for the perturbation 
of m(v.r), 8m(v,r), can be written in the quasi-adiabatic approximation as 



dv 



8m + ATtrX v + 0{Lij), (28) 
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where is the Gaussian stochastic component defined through the noise kernel components i Here one is 

interested in the fluctuations near the apparent horizon. The noise kernel there has not yet been computed, however, 
far from the horizon it has been estimated by Wu and Ford . They found a correlation time for the fluctuations of 
the stress tensor of the order of M and smearing the 2-point functions over this correlation time they found fluctuations 
of the order of 1/M 4 . Form this Hu and Roura [7(3] deduced that 

mtWs-j^Siv-v'), (29) 

which for times larger than the correlation time reproduces the Wu-Ford result (66|. Here £(v) = (47rr 2 ££)[u, r ~ 
6M(v)], relatively far from the horizon. 

One then assumes that as a consequence of the divergenceless property of the stochastic tensor which leads to an 
equation like (|2"5|) , but where replaces (T*), one may be able to connect in a simple way the value of the stochastic 
source near the horizon with its value at large radii. Then, under the assumption that equation Q29p is still valid at 
the horizon Hu and Roura [7(| found that if Mq is the initial mass of the black hole the fluctuations of M(v) are, 

«<-(£)'. W 

which imply that the fluctuations grow in time and they become of the order of the mean value 8M(v) ~ M(v) when 
M(v) ~ Mq . From the evaporation rate equation (|27p the evaporation time of the black hole is of order of Mq and 

2 /3 

thus it reaches the mass M after a significant fraction of the evaporation time. For instance for a black hole with 
an initial mass of the order the solar mass Mq = 1 Mq the fluctuations become comparable to its mean value when 
it reaches a Schwarzschild radius of the order of r$ ~ 10 nm and thus its mass is still much larger than the Planck 
mass and the semiclassical approximation should hold. 

Therefore, according to the previous estimation the fluctuations grow and accumulate in time, the source of this 
accumulation is the non local term in the Einstcin-Langcvin equation, which originates the first term on the right 
hand side of equation (|28[) . This result agrees with the estimation by Bekenstein [67j who also found this long time 
enhancement of the fluctuations. It differs, however, from the estimations by Wu and Ford [66] who neglected the non 
local term. If true, this result seems to point to the breaking of the semiclassical approximation well before the black 
hole reaches the Planck mass and a re-examination of the evaporation process may be needed. At the moment we 
have to take this result with a grain of salt, as some of the approximations used connecting the stochastic tensor near 
the horizon and far from it are not totally correct [70j. But it seems clear that more work is needed, in particular a 
better approximation for the noise kernel near the black hole horizon even within the semi-adiabatic approximation 
is needed. The noise kernel must be treated as a distribution, it is singular in the coincidence limit and for null 
separated points, but it is finite if properly smeared with smooth functions suitably integrated in time as well as in 
space. 



VII. SUMMARY AND OUTLOOK 



We have reviewed the semiclassical theory of gravity as the theory of the interaction of classical gravity with 
quantum matter fields. The most important equations in this theory are the semiclassical Einstein equations |T]) which 
describe the back-reaction of the gravitational fluctuations in its interaction with the quantum fields. We noticed 
that the theory may seriously fail when the fluctuations on the stress tensor of the quantum fields are significant. 
We have then sought an axiomatic approach by which the semiclassical equations can be corrected in order to take 
into account those fluctuations. These equations turn out to be uniquely defined and are the Einstein-Langevin 
equations ((SJ which are linear in the metric perturbations h a b over the semiclassical background. These equations 
predict stochastic fluctuations in the metric perturbations induced by the stress tensor fluctuations described by the 
noise kernel ([3]). We have also noted that the stochastic correlations of the metric perturbations predicted by the 
Einstein-Langevin equations reproduce the quantum metric correlations of the quantum theory of gravity interacting 
with N matter fields, in the large N expansion. 

We have finally used the stochastic theory in the inflationary cosmological context. We have computed the two- 
point correlation functions of the metric fluctuations during a quasi-de Sitter expansion induced by the stress tensor 
fluctuations of the inflaton field. The results are in agreement with other approaches to the same problem [62| . 
an approximate Harrison-Zel'dovich spectrum is predicted. We noticed that in our approach the quantum fields 
and the gravitational fields are treated separately, and this may have some advantages to go one step further and 
consider the quantum field fully, not just to linear order. We have also considered a second application in black hole 
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physics. We have argued that for a large evaporating black hole, fluctuations can accumulate over time and become 
significant before reaching Planck scales. But more work is needed to confirm this calculation and to explore its 
possible consequences. 



Other applications of stochastic semiclassical gravity to semiclassical cosmology have been performed [74| . some 
including thermal fields [H, [7f| . It has been shown that noise produced by a quantum field on the cosmological scale 
factor of an isotropic closed Friedmann-Robertson- Walker, in the presence of a cosmological constant, ma y ta ke the 
scale factor from a region where it is nearly zero to a region where it describes a de Sitter inflationary era [761 ] . Thus 
jumping over the barrier by activation, this is the semiclassical analogue of the tunneling from nothing in quantum 
cosmology [zj [zS [zH and gives yet another mechanism to produce inflation. Finally, stochastic gravity has also been 
used to formulate a criteria for the validity and stability of semiclassical gravity [I?} . In particular it has been shown, 
in this context, that flat spacetime as a background solution of semiclassical gravity is stable. 
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